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The  Army  High  Performance  Comput¬ 
ing  Research  Center,  a  collaboration 
between  the  U.S.  Army  and  a  consor¬ 
tium  of  university  and  industry  partners, 
develops  and  applies  high  performance 
computing  capabilities  to  address  the 
Army's  most  difficult  scientific  and  engi¬ 
neering  challenges. 

AHPCRC  also  fosters  the  education  of 
the  next  generation  of  scientists  and 
engineers— including  those  from  racially 
and  economically  disadvantaged  back¬ 
grounds— in  the  fundamental  theories 
and  best  practices  of  simulation-based 
engineering  sciences  and  high  perfor¬ 
mance  computing. 

AHPCRC  consortium  members  are: 
Stanford  University,  High  Performance 
Technologies  Inc.,  Morgan  State  Univer¬ 
sity,  New  Mexico  State  University  at  Las 
Cruces,  the  University  of  Texas  at  El  Paso, 
and  the  NASA  Ames  Research  Center. 


Strong,  impact-resistant  materials 
lighten  the  soldiers  load,  give  the 
soldier  increased  protection,  and 
minimize  risk  of  injury.  Computer 
simulation  allows  designers  to 
compare  and  evaluate  numerous 
mechanical  and  material 
configurations  without  the  time 
and  expense  of  building  physical 
models.  Complex,  realistic  computer 
simulations  enable  more  sophisticated 
and  intentional  approaches  to  material 
design  and  better  understanding  of  a 
material’s  modes  of  deformation  and 
failure. 


Computer  simulations, 
clockwise  from  top  left:  Ballistic 
impact  on  protective  fabric; 
ballistic  impact  on  a  soft 
target;  air  flow  in  the  wake  of  a 
flapping  wing;  airfoil  pressure 
contour  map.  (All  graphics 
courtesy  of  the  researchers.) 


Soldier  survivability  is  also  increased 

by  using  mechanical  devices,  such  as  automated  aerial  drones,  to 
stand  in  for  humans  where  hazardous  or  tedious  work  is  involved. 
Computer  simulations  assist  in  developing  mechanical  designs  that 
exploit  solutions  already  in  use  by  biological  organisms.  Streamlined, 
massively  parallel  high  performance  computing  structural  codes 
allow  researchers  to  examine  many  relevant  physical  factors 
simultaneously,  or  to  examine  parameters  in  isolation  (often  not 
possible  in  the  physical  world),  thus  creating  models  that  more 
closely  simulate — or  improve  on —  real-world  phenomena. 


In  this  issue  of  the  AHPCRC  Bulletin,  we  introduce  the  projects  in 
AHPCRC  s  Lightweight  Combat  Systems  Survivability  technical  area 
and  the  researchers  working  on  these  projects.  Also  inside:  awards 
and  recognition  received  by  AHPCRC  consortium  partners,  and  a 
list  of  publications  and  presentations  based  on  AHPCRC  research.  A 
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Multiscale  Ballistic  Fabric  Nets 

Body  armor  ranks  with  water,  ammunition,  and 
weapons  as  one  of  the  heaviest  items  worn 
or  carried  by  troops,  according  to  engineers 
on  the  Ballistics  Technology  Team  at  the  U.S.  Army 
Soldier  Systems  Center  (Natick,  MA).  Protecting 
soldiers  and  their  vehicles  and  equipment  from  flying 
projectiles  requires  shielding  that  resists  penetration 
and  that  is  lightweight  enough  to  not  interfere  with 
normal  activities  and  operations.  In  addition,  shield¬ 
ing  material  must  retain  its  effectiveness  even  after 
long  periods  in  storage  and  exposure  to  light,  heat, 
and  humidity. 

High  performance  computing  simulations  allow  re¬ 
searchers  to  develop  better  ballistic  shielding  materials 
by  examining  the  effects  of  ballistic  impact  frame-by- 
frame.  Researchers  can  use  simulations  to  test  various 
armor  components  singly  or  in  groups.  Simulations 
may  suggest  effective  material  configurations  that  are 
not  intuitively  obvious  from  experimental  data  alone. 
After  simulations  narrow  the  field  to  the  most  promis¬ 
ing  armor  configurations,  they  can  be  fabricated  and 
tested  under  real-life  conditions. 

AHPCRC  primary  investigators  Tarek  Zohdi  (associ¬ 
ate  professor,  mechanical  engineering,  University  of 
California,  Berkeley)  and  Charbel  Farhat  (professor, 
aeronautics,  astronautics,  and  mechanical  engineer¬ 
ing,  Stanford  University),  research  associate  Philip 
Avery,  and  postdoc  David  Powell  (both  of  Stanford) 
are  bringing  ballistic  fabric  modeling  into  the  high 
performance  computing  arena  by  adapting  existing 

Time-lapse  computer  simulation  of  fabric 
deformation  and  rupture. 

I 


codes  for  simulation  and 
finite  element  analysis  to 
run  in  parallel  processing 
environments. 

Because  there  has  not  been 
an  abundance  of  new  ma¬ 
terials  to  test,  recent  efforts 


Ballistic  impact  produces  multiple  physical  effects  at 
high  speeds — understanding  these  effects  will  assist  in 
designing  protective  fabrics  that  can  withstand  im¬ 
pact  by  various  types  of  projectiles.  Using  computer 
simulations  to  model  the  behavior  of  a  fabric  under 
impact  from  a  high-speed  bullet  or  shrapnel  fragment 
requires  the  capability  to  model  the  nonlinear  solid 
dynamics  typical  of  materials  undergoing  rapid  defor¬ 
mation  caused  by  localized  impact. 

In  addition  to  adapting  existing  software  codes,  Zohdi, 
Farhat,  Avery,  and  Powell  are  building  new  capabili¬ 
ties  into  their  modeling  and  simulation  software  to 
account  for  such  rapid  deformation  effects.  They  are 
adding  simulations  for  various  projectile  shapes,  ac¬ 
counting  for  imperfections  incorporated  during  the 
weaving  of  protective  fabrics,  examining  the  way  the 
protective  fabric  is  attached  to  the  underlying  struc- 


have  focused  on  optimizing 
the  performance  of  existing 
materials.  Fiber  strength  is 
important,  but  fabric  geom¬ 
etry  also  matters,  because 
projectiles  can  penetrate  by 
pushing  fibers  aside,  even  if 
the  fibers  remain  unbroken. 

Factors  that  affect  perfor¬ 
mance  of  ballistic  fabrics 
include  the  thickness  and 
strength  of  the  fibers,  weav¬ 
ing  patterns,  and  the  method  of  attaching  the  fabrics 
to  their  substrates  (pinned  at  the  corners  vs.  fastened 
along  the  sides,  for  example).  Customarily,  the  Army 
has  evaluated  ballistic  fabrics  by  making  and  testing 
physical  prototypes.  Testing  every  combination  of 
characteristics  to  find  optimal  combinations  is  time- 
consuming  and  expensive,  however. 


Zylon  fabric  test  samples  are 
ruptured  using  a  cylindrical 
penetrator.  (Photo  and 
graphics  courtesy  David 
Powell,  Stanford.) 
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ture,  and  simulating  the  propagation  and  growth 
of  flaws  in  fabrics.  They  have  begun  preliminary 
work  on  modeling  fiber-based  composite  mate¬ 
rials  in  order  to  explore  the  properties  of  these 
complex  systems. 


Putting  the  Model  to  the  Test 

The  simulation  method  developed  for  this 
project  represents  ballistic  fabrics  from  the  fibril 
level  up,  so  that  the  same  method  can  be  applied 
to  many  different  materials,  weaving  types,  and 
modes  of  attachment  to  a  substrate.  For  both  the  simu¬ 
lation  and  the  real-world  fabrics,  thin  filaments  called 
fibrils  are  twisted  together  into  strands  of  yarn,  which 
are  then  woven  into  fabric  sheets.  Fibril  properties  are 
relatively  simple  to  model,  they  can  be  obtained  readily 
from  the  manufacturer,  and  they  are  simple  to  test  and 
verify. 

The  model  addresses  tensile  strains  only,  as  tension 
(rather  than  compression)  is  most  important  when 
assessing  rupture  properties  due  to  ballistic  impact. 
Because  a  fibril  is  very  thin  relative  to  its  length,  tensile 
deformation  can  be  described  as  a  one- dimensional 
response  to  an  applied  uniaxial  stress.  Axial  strain  for 
the  fabric  is  assumed  to  lie  in  the  region  between  2% 
and  10%  prior  to  rupturing.  Zylon  fabric,  which  was 
used  in  the  experiments  with  which  the  simulations 
were  compared,  ruptures  at  3%  strain.  Zylon  yarn  has 
350  fibrils  per  strand,  and  the  simulated  yarn  response 
is  obtained  by  summing  the  responses  for  all  the  fibrils. 

Misalignment,  in  which  the  directions  of  the  fibril  axes 
exhibit  a  statistical  distribution,  occurs  during  manu¬ 
facturing  (indeed,  it  is  almost  impossible  to  avoid). 

This  misalignment  is  a  beneficial  property  of  the  yarn. 
Individual  fibrils  fail  at  different  times  in  response  to 
stress,  because  of  the  variation  in  their  orientation  with 
respect  to  that  stress.  This  causes  a  given  yarn  strand  to 
fail  gradually  rather  than  catastrophically.  Differences 
in  yarn  strands  tend  to  be  small,  because  the  large 
number  of  fibrils  in  each  strand  average  out  statisti¬ 
cally,  but  the  misalignment  effect  is  a  necessary  com¬ 
ponent  of  a  good  macroscale  model. 


l  i 

Left:  Conventional  contact  model  performs  calculations  for 

each  body  using  separate  processors.  Right:  Areas  in  close 
physical  contact  are  processed  together,  reducing  data 
transfer  between  processors. 

work  of  nonlinear  trusses  (yarn)  pinned  together  at 
nodes  where  two  trusses  meet.  (“Nonlinear”  refers 
to  materials  that  do  not  experience  strain  in  propor¬ 
tion  to  the  amount  of  stress  applied.)  This  model  was 
constructed  using  FEM,  a  finite  element  modeling 
program  used  at  Stanford.  Sandia  Laboratories’  ACME 
contact  library  was  added  to  FEM  to  supply  appropri¬ 
ate  algorithms  to  search  for  contact  entities  and  to 
enforce  contact  constraints. 

When  simulations  deal  with  parts  of  various  bodies 
(e.g.,  a  projectile  or  a  sheet  of  fabric)  that  are  in  physi¬ 
cal  contact,  such  as  during  an  impact  event,  significant 
amounts  of  information  are  passed  back  and  forth  be¬ 
tween  computer  processors,  which  increases  the  time 
needed  to  complete  the  calculations.  Optimizing  the 
performance  of  the  contact  search  proved  to  be  depen¬ 
dent  on  the  choice  of  subdomains,  because  separate 
processors  perform  calculations  for  each  subdomain. 
Typically,  subdomains  group  together  elements  that 
are  in  the  same  body.  A  more  efficient  method  defines 
subdomains  as  those  regions  which  are  in  close  spatial 
proximity  to  each  other,  regardless  of  which  body  they 
are  in  (see  figure  above). 

The  problem  of  a  small  projectile  striking  a  single  fab¬ 
ric  sheet  presents  an  area  of  contact  that  involves  only 
a  few  elements.  This  makes  it  difficult  to  divide  the 
workload  over  more  than  about  10  processors.  At  the 
other  extreme,  when  two  sheets  come  together  and  the 


The  AHPCRC  model  treats  a  woven  fabric  as  a  net- 
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Simulating  Fracture  and 
Penetration 

nderstanding  the  effects  of  ballistic  impact 
of  ammunition  or  shrapnel  on  human  soft 
tissues  has  the  potential  to  aid  medical  pro¬ 
fessionals  in  treating  battlefield  injuries,  planning  re¬ 
constructive  surgery,  and  understanding  more  about 
wound  healing  processes.  In  addition,  this  knowledge 
can  be  used  to  evaluate  various  materials  and  designs 
for  ammunition  and  armor.  Human  testing  of  this 
type  is  of  course  impractical.  Blocks  of  alternative  tis¬ 
sue  simulants,  called  ballistic  gels,  are  used  instead. 

Constructing  computer  simulations  of  material 
behavior  under  shock  or  ballistic  impact  is  especially 
difficult  because  complex  changes  occur  rapidly.  A 
soft  material  such  as  a  polymer  or  ballistic  gel  absorbs 
and  dissipates  energy  from  the  impact,  deforms, 
and  cracks,  in  less  than  a  few  milliseconds.  Build¬ 
ing  a  computer  simulation  that  captures  all  of  these 
changes  requires  repeatedly  solving  millions  of  equa¬ 
tions  to  create  a  series  of  “snapshots”  taken  at  the  rate 
of  millions  of  frames  per  second  and  extending  over 
several  milliseconds.  Each  change  affects  the  changes 
after  it,  and  physical  effects  interact  to  alter  the  course 
of  events.  In  addition,  each  snapshot  comprises  mil¬ 
lions  of  small  sections  (finite  elements)  of  a  structure, 
each  interacting  with  and  affecting  the  behavior  of 
the  adjacent  sections.  The  number  of  computations 
required  to  construct  a  realistic  simulation  can  tax  the 
resources  of  even  the  best  high  performance  comput¬ 
ers.  Thus,  it  is  necessary  to  frame  the  problem  and 
program  the  computational  codes  in  a  way  that  uses 
computing  resources  efficiently,  without  sacrificing 
accurate  results. 

AHPCRC  researcher  and  assistant  professor  Adrian 
Lew  and  graduate  students  Raymond  Ryckman  and 
Ramsharan  Rangarajan  (mechanical  engineering 
department,  Stanford  University),  are  working  with 
AHPCRC  staff  scientist  Mark  Potts  (High  Perfor¬ 
mance  Technologies,  Inc.)  to  create  and  adapt  meth¬ 
ods  for  modeling  the  evolution  of  domains  (regions 


Simulation  of  a  spherical  projectile  penetrating  a  ballistic  gel. 
(Photo  and  graphics  courtesy  Adrian  Lew,  Stanford.) 


with  uniform  properties)  and  crack  paths  so  that  they 
function  in  a  parallel  processing  environment. 

The  group  is  devising  methods  to  simulate  soft  ma¬ 
terials  undergoing  fracture  and  penetration.  When 
they  began  their  work,  there  were  no  existing  simula¬ 
tion  methods  that  were  sufficiently  robust  to  carry 
out  accurate  calculations  over  the  time  scale  of  an 
actual  impact  event.  Because  of  the  computational 
resources  needed  for  such  an  endeavor,  any  computer 
codes  would  have  to  be  scalable;  that  is,  the  codes 
would  need  to  run  well  on  any  number  of  available 
processors.  The  group  needed  to  develop  algorithms 
for  complex  motions  and  interactions  between  the 
projectile  and  the  target  that  would  work  well  in  such 
an  environment. 

At  present,  a  model  is  being  constructed  for  the 
effects  of  a  bullet  hitting  a  ballistic  gel.  Ballistic  gel 
simulations  under  development  include  shock, 
penetration,  and  contact  effects,  and  the  way  the  gel 
absorbs  energy  from  an  impact.  Lew  foresees  future 
efforts  that  incorporate  aspects  of  material  modeling: 
how  does  impact  with  tissue  or  a  polymer  differ  from 
impact  with  wood,  metal,  or  concrete?  What  happens 
when  friction-generated  heat  partially  melts  the  poly¬ 
mer  or  gel?  What  are  the  effects  of  a  projectile  tum¬ 
bling  within  a  fissure?  How  does  a  network  of  cracks 
form  and  spread?  The  Stanford  group  is  planning  to 
compare  their  simulations  with  experimental  studies 
conducted  at  the  University  of  California  at  Berkeley 
at  a  later  stage  of  the  project. 
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Breaking  the  Lock-Step 

The  AHPCRC  group  is  currently  working  on  a  par¬ 
allel  code  for  nonlinear  solid  dynamics  that  incorpo¬ 
rates  an  algorithm  capable  of  asynchronously  con¬ 
structing  the  part  of  each  snapshot  corresponding 
to  different  spatial  locations  (parallel  asynchronous 
time  integration,  or  pAVI).  This  approach  allows 
each  section  of  a  large  problem  to  be  computed  at 
the  speed  and  degree  of  detail  that  best  suits  the 
complexity  of  that  section.  This  differs  from  con¬ 
ventional  problem-solving  approaches,  in  which 
differential  equations  for  each  of  the  millions  of 
elements  in  a  problem  are  integrated  using  the  same 
time  increments,  using  the  same  number  of  proces¬ 
sors.  Potts  explains:  “With  the  lock  step  approach, 
lots  of  computational  time  ends  up  getting  devoted 
to  solving  nondescript  sections  of  the  problem  using 
very  small  time  steps  where  they  otherwise  could  be 
solved  using  much  larger  time  steps.” 

For  problems  in  which  some  sections  are  especially 
complex  (e.g.,  complicated  geometries,  rapid  shape 
changes),  the  time  increments  between  individual 
calculations  must  be  very  small  in  order  to  capture 
the  requisite  degree  of  detail — in  much  the  same  way 
as  high-speed  video  uses  thousands  of  frames  per 
second  to  capture  the  flight  of  a  bullet  or  the  impact 
of  a  falling  drop  of  water.  Dividing  the  complex  parts 
of  the  problem  among  several  processors  allows  for 
frequent  updates  without  slowing  down  the  overall 
process.  Simpler  or  less  important  sections  of  the 
problem  that  do  not  require  this  degree  of  detail  can 
proceed  with  fewer  updates  (larger  time  increments 
between  calculations)  and  occupy  fewer  processors. 


Photograph  of  a  projectile  path  through  a  ballistic  gel. 


Thus,  the  pAVI  approach 
reduces  the  need  for 
trade-offs  between  ac¬ 
curacy  and  efficiency. 
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Redrawing  the 
Boundaries 

The  mathematical  ap¬ 
proach  taken  by  the 
AHPCRC  group  is  based  Top:  A  boundary-fitting  mesh, 
on  numerical  methods  Bottom:  A  mesh  suitable  for 

for  solving  partial  differ-  use  with  immersed  boundary 
,  methods, 

ential  equations,  which 

track  physical  quantities 

that  change  continuously  over  time.  A  mathematical 
representation  of  a  soft  structure,  such  as  a  ballistic 
gel,  is  divided  into  many  small  elements,  each  of 
which  changes  shape  in  response  to  the  simulated 
ballistic  impact.  As  these  elements  stretch  and  bend 
away  from  their  initial  shape,  the  mathematics 
describing  them  becomes  more  complex.  At  some 
point,  the  calculations  must  pause,  and  a  process 
known  as  “remeshing”  draws  new  reference  gridlines 
onto  the  deformed  soft  body,  based  on  localized 
displacement  and  velocity  fields.  This  procedure  may 
need  to  be  repeated  many  times  during  the  course 
of  a  simulation,  and  it  slows  the  calculation  process 
considerably. 


The  Stanford  group  uses  an  adaptive  remeshing 
algorithm  based  on  an  immersed  boundary  method 
(see  figure  above).  This  method,  which  is  often  used 
to  simulate  solid  bodies  moving  through  fluids  (or 
gels,  in  this  case),  uses  a  mesh  that  encompasses  the 
domain  of  the  problem  while  only  approximately 
following  the  boundaries  between  the  soft  body  and 
the  solid  body  passing  through  it.  This  makes  gener¬ 
ating  the  mesh  much  easier  at  the  expense  of  com¬ 
plicating  the  imposition  of  boundary  conditions. 
Implementing  this  concept  for  a  three-dimensional 
problem  is  one  of  the  major  areas  of  research  for 
this  project,  because  of  the  ability  of  this  method  to 
overcome  many  of  the  difficulties  associated  with 
automatic  mesh  generation. 
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Optimizing  Lift  vs.  Drag  for  Small 
Airfoils 


Eagles  soar,  hummingbirds  whir.  How  a  bird 
flies  depends  largely  on  how  big  the  bird  is. 
Likewise,  hummingbird- sized  aerial  surveil¬ 
lance  and  reconnaissance  drones  cannot  depend  on 
the  same  wing  configurations  and  aerodynamics 
as  their  eagle-sized  counterparts.  As  aerial  vehicles 
become  progressively  smaller,  the  viscosity  of  the  air 
takes  on  a  greater  importance  ( see  “All  About  Eddy  ” 
p.  7).  The  ratio  of  the  vehicles  inertial  force  (propel¬ 
ling  it  forward)  to  the  viscous  force  of  the  air  (hold¬ 
ing  it  back)  decreases.  To  keep  a  tiny  vehicle  airborne 
on  just  the  amount  of  fuel  that  it  can  carry  requires 
wings  that  can  generate  the  most  lift  with  the  least 
drag. 


AHPCRC  researchers  Professor  Antony  Jameson 
(Stanford  University,  Aeronautics  and  Astronautics) 
and  his  student  Matt  Culbreth  are  using  massively 
parallel  high  performance  computing  simulations  to 
optimize  airfoil  shapes  and  to  identify  the  attainable 
performance  limits  for  a  given  vehicle  size  and  con¬ 
figuration.  For  a  given  amount  of  lift,  the  computa¬ 
tions  test  various  shapes  to  see  which  shapes  produce 
the  least  amount  of  drag.  By  using  computational 
models,  the  researchers  can  test  a  large  number  of 
vehicle  designs  much  more  rapidly  than  if  they  had 
to  generate  physical  prototypes  for  each  design. 


This  page:  Two  examples  of  airfoils  generated  by 
Jameson's  group  at  Stanford.  Colored  areas  represent 

pressure  contours. 

The  Stanford  researchers  are  performing  their  cal¬ 
culations  using  low  values  (1000  to  10,000)  for  the 
ratio  of  inertial  to  viscous  forces,  called  the  Reyn¬ 
olds  number.  Culbreth  and  Jameson  have  generated 
several  airfoil  configurations  using  a  small  number 
of  geometry  control  points.  Their  work  has  produced 
airfoil  geometries  corresponding  to  a  given  pressure 
distribution,  and  they  have  successfully  found  the 
minimum  amount  of  drag  for  a  fixed  lift  coefficient. 

The  AERO-F  simulation  code  that  they  use  is  mas¬ 
sively  parallel;  that  is,  it  runs  on  computer  systems  in 
which  the  work  is  distributed  across  many  individual 
nodes.  Each  node  consists  of  at  least  one  proces¬ 
sor,  its  own  memory,  and  a  link  to  the  network  that 
connects  all  the  nodes  together.  The  largest  super¬ 
computers  have  several  hundreds  of  thousands  of 
nodes  in  the  same  “box,”  but  parallel  systems  can  also 
consist  of  clusters  of  conventional  computers  linked 
in  such  a  way  that  they  pass  information  back  and 
forth  as  seamlessly  as  a  single  supercomputer.  AERO- 
F  is  also  scalable,  making  it  possible  to  run  equally 
well  on  systems  of  various  sizes. 

Computer  simulations  have  already  identified  prom¬ 
ising  airfoil  shapes  that  run  counter  to  an  engineers 
normal  intuition — and  that  Nature  never  thought  of. 
Camber  (the  difference  in  curvature  between  the  top 
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and  bottom  surfaces  of  an  airfoil)  can  be  optimized 
to  reduce  drag  or  to  increase  the  angle  at  which  the 
aircraft  begins  to  stall.  If  the  curvature  is  greater  on 
the  top  surface,  the  camber  is  positive,  while  a  greater 
curvature  on  the  bottom  surface  produces  negative 
camber — a  configuration  known  as  a  supercritical 
airfoil  shape,  which  has  been  used  to  improve  the  lift- 
to-drag  ratio  for  high-speed  aircraft. 

Micro  aerial  vehicles  are  typically  built  for  minia¬ 
turization  rather  than  speed,  and  flapping  wings  are 
a  necessity  to  generate  enough  lift  to  remain  aloft. 
Simulations  have  identified  combinations  of  flap¬ 
ping  frequency  and  amplitude  that  do — and  don’t — 
generate  sufficient  lift.  The  simulations  show  a  trail  of 
vortices  behind  the  flapping  wing,  which  can  either 
help  the  vehicle  (or  bird)  stay  aloft  or  hinder  it.  ( See 
next  article  for  more  about  this  work.) 

Recent  work  is  moving  toward  a  3D  unsteady  optimi¬ 
zation  of  flapping  motions  for  hovering  and  forward 
flight.  The  goal  is  to  couple  an  advanced  computation¬ 
al  fluid  dynamics  code  such  as  AERO-F  with  the  op¬ 
timization  libraries  developed  elsewhere  in  AHPCRC 
in  order  to  synthesize  wing  motions  and  deformations 
that  optimize  a  specific  performance  metric. 

Several  iterations  of  a  base  mesh  have  been  generated 
for  the  initial  optimizations,  incorporating  refine¬ 
ments  to  improve  accuracy,  convergence,  and  resolu¬ 
tion  of  flow  features,  while  attempting  to  minimize 
the  computational  cost  as  much  as  possible.  ( See 
previous  article  for  more  about  meshing  schemes.) 

The  appropriate  solver  parameters  are  being  de¬ 
termined  for  the  AERO-F  code  running  unsteady 
simulations  at  Reynolds  numbers  between  1,000  and 
10,000.  This  process  involves  determining  how  to 
automate  flow  solutions  and  the  post-processing  of 
results  so  that  they  can  be  integrated  with  the  optimi¬ 
zation  codes.  At  present,  simulations  involving  wing 
pitching  and  plunging  are  being  prepared  to  assess 
the  trade-offs  between  mesh  density  and  the  ability  to 
capture  vorticity.  ★ 


All  about  Eddy 

Airfoil  modeling  is  not  a  new  science,  although 
steady  advancements  in  computing  power  enable 
the  creation  of  simulations  with  an  unprecedent¬ 
ed  degree  of  complexity.  Several  methods  exist 
for  expressing  fluid  flow  in  mathematical  terms. 

Large  eddy  simulation  (LES)  is  a  technique  for  solv¬ 
ing  the  differential  equations  associated  with  tur¬ 
bulent  flow.  The  system  of  airfoil  and  flowing  fluid 
is  represented  by  points  on  a  grid,  and  calculations 
are  carried  out  at  these  points.  A  finer-meshed 
grid  produces  a  greater  accuracy,  but  is  more 
expensive  in  terms  of  computer  resources  and 
time  requirements.  In  the  LES  method,  simula¬ 
tions  are  generated  for  only  the  larger  eddies,  and 
the  smaller  eddies  are  handled  using  a  sub-grid 
scale  model.  LES  represents  an  intermediate  level 
between  direct  numerical  simulation  (DNS),  which 
models  turbulent  motion  at  all  relevant  scales, 
and  the  more  approximate  Reynolds-averaged 
Navier-Stokes  (RANS)  model. 

In  variational  multiscale  LES,  calculations  are 
separated  into  two  or  more  groups  that  address 
different  size  scales.  The  filters  used  in  the  LES 
method  are  replaced  by  variational  projections, 
an  approach  that  produces  greater  efficiency  and 
accuracy  than  the  conventional  LES  method. 

Detached  eddy  simulations  (DES)  uses  a  RANS 
method  to  solve  regions  that  have  less  than  a 
specified  turbulent  length  scale  and  that  are  near 
solid  boundaries.  If  the  turbulent  length  scale 
exceeds  the  grid  dimension,  the  regions  are  solved 
using  LES.  The  DES  method  requires  less  grid  reso¬ 
lution  than  the  LES  method,  which  reduces  the 
cost  of  the  computation. 

Source:  Adapted  from  Wikipedia. 
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Flapping  and  Twisting  Wings 

Building  a  robotic  hummingbird  is  probably 
not  the  first  thing  one  associates  with  Army 
research,  but  bird-sized  robotic  aerial  vehicles 
are  much  in  demand.  Micro  aerial  vehicles,  or  MAVs, 
could  perform  surveillance,  hazardous  materials 
detection,  terrain  mapping,  or  data  relay  functions  in 
situations  where  it  would  be  hazardous  or  impracti¬ 
cal  to  send  humans.  MAVs  don’t  need  sleep,  and  they 
don’t  get  bored  by  looking  at  the  same  terrain  for 
hours  at  a  stretch.  Unmanned  aerial  vehicles  (UAVs), 
larger  counterparts  to  the  MAVs,  have  proven  useful 
for  military  operations  in  Afghanistan,  but  the  current 
generation  of  vehicles  requires  a  runway  for  takeoff, 
and  they  are  unable  to  hover. 

AHPCRC  researchers  from  New  Mexico  State  Uni¬ 
versity  are  working  on  a  robotic  hummingbird  wing 
to  gain  insight  into  how  to  mimic  a  hummingbird’s 
maneuverability.  They  have  incorporated  information 
from  videos,  provided  by  Bret  Tobalske  (then  at  the 
University  of  Portland,  now  at  the  University  of  Mon¬ 
tana),  of  the  vortices  that  live  hummingbirds  make  as 
they  fly.  Hummingbird  wings  flap  and  twist,  creating 
eddy  currents  that  give  the  bird  lift  and  momentum. 

Mingjun  Wei  and  Banavara  Shashikanth,  NMSU  assis¬ 
tant  professors  of  mechanical  and  aerospace  engineer¬ 
ing,  head  up  the  project.  The  late  James  Allen,  NMSU 
assistant  professor  of  mechanical  engineering,  over¬ 
saw  the  project  in  its  early  stages.  Prof.  Tom  Burton, 
department  head,  was  instrumental  in  bringing  the 
NMSU  group  into  the  AHPCRC  program,  and  he  par¬ 
ticipates  in  the  research  group  meetings.  Prof.  Charbel 
Farhat  of  Stanford  University  is  providing  additional 
computer  simulation  work  on  lift  and  drag.  NMSU 
graduate  students  Tao  Yang,  Humberto  Bocanegra 

Evans,  Mohammad 
Mazharul  Islam, 
and  Ramiro  Chavez 
receive  AHPCRC 
support. 

3:1  scale  model  of  a 
hummingbird  wing. 


The  research  team  is 
building  computer 
simulations  that  will 
help  them  understand 
and  imitate  the  com¬ 
plex  wing  motions 
that  come  naturally  for 
small  birds  and  in¬ 
sects.  At  the  scale  of  a 
small  bird,  you  have  to 
deal  with  factors  that 
don’t  even  register  for 
large  aircraft  (or  even 
large  birds) — a  gust 
of  air  from  a  building’s 
exhaust  vent,  or  drag 
generated  by  the  viscos¬ 
ity  of  the  air,  for  ex¬ 
ample.  Birds  and  insects  have  been  dealing  with  wind 
gusts,  fuel  intake  and  consumption,  lift  and  drag,  and 
obstacle  avoidance  issues  for  millennia.  They  fly  under 
all  sorts  of  weather  conditions,  alone  or  in  flocks  or 
swarms.  (See  previous  article  for  more  on  small  bird 
flight.) 


Mechanical  wing  testing  in  a 
low-turbulence  flow  tank  at 
NMSU.  (Photos  and  graphics 
courtesy  Mingjun  Wei,  Banavara 
Shashikanth,  NMSU.) 


How  does  a  bird  do  it? 

Among  the  problems  an  MAV  designer  has  to  face 
is  how  to  gather  and  store  enough  energy  to  sustain 
flight  without  making  the  vehicle  too  heavy  to  fly.  The 
hummingbird  does  this  by  feeding  every  10-15  min¬ 
utes  during  its  waking  hours  and  slipping  into  a  deep 
sleep  state  at  night.  It  would  be  impractical  to  operate 
an  MAV  on  this  schedule,  so  efficient  energy  usage  is 
key  to  the  success  of  these  vehicles. 


Birds  and  insects  make  very  little  noise  as  they  fly, 
and  they  recover  quickly  if  they  graze  a  tree  branch  or 
other  obstacle.  These  will  be  important  features  in  de¬ 
signing  MAVs,  along  with  agility  and  maneuverability. 
Hummingbirds  are  especially  agile — they  can  take  off 
from  a  stationary  position,  hover,  and  change  direc¬ 
tions  in  midair.  During  flight,  they  sweep  their  wings 
forwards  and  backwards  while  plunging  (oscillating) 
and  pitching  (rotating)  their  wings,  tracing  figure- 
eight  paths  with  their  wing  tips. 
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Building  a  better  bird 

The  earthbound  engineers  at  NMSU  are  using  high- 
performance  computer  (HPC)  simulations  to  explore 
the  best  combinations  of  plunging,  pitching,  and 
twisting  motions  for  producing  thrust.  They  are  test¬ 
ing  various  frequencies  and  amplitudes  of  the  flapping 
and  twisting  motions  to  determine  what  works  best. 
They  are  developing  a  three-dimensional  HPC  version 
of  a  code  based  on  a  modified  immersed  boundary 
method  (see  page  5),  which  models  a  flowing  fluid 
interacting  with  a  flexible,  elastic  structure  (such  as  a 
wing)  by  computing  the  flow  and  the  structure  simul¬ 
taneously  on  the  same  mesh  grid  in  a  fully- coupled 
fashion.  Dynamical  models  for  flow-structure  interac¬ 
tions  are  being  developed  for  some  simple  configura¬ 
tions.  Wei  and  Yang  are  doing  numerical  simulations, 
and  Shashikanth  and  Islam  are  developing  theoretical 
models. 

Evans,  Chavez,  and  their  co-workers  have  built  a 
model  wing  that  they  are  testing  in  a  low- turbulence 
water  channel,  and  they  will  use  their  findings  to 
validate  and  improve  their  computer  models.  This  will 
enable  them  to  find  useful  configurations,  maybe  even 
ones  that  Nature  hasn’t  thought  of  yet,  and  save  their 
mechanical  model-building  resources  for  the  best  can¬ 
didates.  Preliminary  runs  have  been  conducted  on  a 
mechanical  model  wing  with  two  degrees  of  freedom 
of  motion,  and  the  concept  design  has  been  completed 
for  a  wing  with  three  degrees  of  freedom. 

Farhat  s  group  at  Stanford  is  pursuing  an  alterna¬ 
tive  simulation  approach  based  on  a  combination  of 
higher-order  ALE  (Arbitrary  Lagrangian  Eulerian) 
and  embedded  methods.  This  approach  is  especially 
useful  for  modeling  the  interactions  of  moving  flexible 
wings  with  air  streams,  because  it  enables  calculations 
based  on  the  interaction  of  both  solid  and  fluid  refer¬ 
ence  frames. 

To  this  effect,  they  have  developed  an  upgraded 
version  of  the  massively  parallel  AERO-F  computer 
code  that  features  a  number  of  explicit  and  implicit 
schemes  that  satisfy  their  discrete  geometric  conser¬ 
vation  laws  (DGCLs) — a  set  of  conditions  imposed 
on  a  moving  computational  grid  in  order  to  preserve 


the  accuracy  and  stability  of  the  calculations  when 
integrated  over  a  period  of  time.  The  formal  accuracy 
of  these  schemes  was  proven  to  be  identical  to  their 
counterparts  on  fixed  grids.  The  team  has  analyzed  the 
amount  of  thrust  generated  by  high-frequency  plung¬ 
ing  motions  using  AERO-F  s  large  eddy  simulation 
(LES)  capabilities.  More  recently,  they  simulated  flap¬ 
ping  motions  (active  plunging  coupled  with  passive 
pitching)  for  various  frequencies. 

In  particular,  LES  runs  were  performed  at  Stanford  on 
massively  parallel  processors  to  investigate  numeri¬ 
cally  the  existence  of  a  critical  frequency,  below  which 
drag  is  generated  and  above  which  thrust  is  produced. 
The  simulations  clearly  demonstrated  the  existence 
of  such  a  critical  frequency.  Previous  experimental 
studies  had  shown  a  surprising  onset  of  nonsymmet- 
ric  wake  patterns  at  higher  frequencies.  The  Stanford 
simulations  also  predicted  these  patterns,  which  pro¬ 
duce  both  lift  and  thrust. 

Insights  gained  through  computational  modeling  and 
simulation  provide  valuable  guidance  in  choosing 
configurations  and  operational  parameters  for  real- 
world  mechanical  wings,  and  for  understanding  the 
factors  that  go  into  making 
a  wing  that  works.  Back  at 
NMSU,  work  continues  on 
improving  the  physical  wing 
models  for  testing  in  the 
water  channel.  The  research 
group  also  plans  to  build 
a  computational  model 
of  a  micro  aerial  vehicle 
with  flapping  wings.  They 
will  test  their  models  us¬ 
ing  plunging,  pitching,  and 
twisting  motions  alone  and 
in  combination  to  deter¬ 
mine  which  motions  pro¬ 
duce  the  most  thrust.  A 


Vortices  created  by  a  flapping 
flexible  airfoil  in  an  incoming 
flow,  modeled  over  a  period  of 
0.5  second. 
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AHPCRC  Tech  Area  1  Research  Team 

Tech  Area  Lead:  Charbel  Farhat 

Principal  Investigators: 


Multiscale  Ballistic  Fabric  Nets 

Tarek  I.  Zohdi 
Mechanical  Engineering 
Department 

University  of  California  Berkeley 

zohdi@newton.berkeley.edu 

(510)-642-9172 


Multidisciplinary  Parametric 
Modeling  and  L/D 
Quantification  and  Optimization 

Antony  Jameson 
Aeronautics  and  Astronautics 
Department 
Stanford  University 
jameson@baboon.stanford.edu 
(650)  725-6208 


and 

Charbel  Farhat 
Mechanical  Engineering 
JP  Department 

Stanford  University 
Ik  cfarhat@stanford.edu 
(650) 723-3840 


Flapping  and  Twisting  Aeroelastic 
Wings  for  Propulsion 

Mingjun  Wei 

Mechanical  and  Aerospace 
Engineering  Department 
New  Mexico  State  University 
mjwei@nmsu.edu 
(575)  646-6034 


Simulation  of  Fracture  and 
Penetration 

Adrian  Lew 

Mechanical  Engineering 
Department 
Stanford  University 
I  e  wa  @  sta  nf  o  rd .  ed  u 
(650)  725-3585 


I  and 

y  ®  Banavara  N.  Shashikanth 

Mechanical  and  Aerospace 
Engineering  Department 
New  Mexico  State  University 
i||^L  shashi@nmsu.edu 

(575)  646-4348 
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AHPCRC  Staff  News 


Dick  Pritchard, 

AHPCRC  Program 
Manager 

Since  2008,  Dr.  Richard 

(Dick)  Pritchard  has  been  Program  Manager 
for  AHPCRC.  He  is  responsible  for  the 
infrastructure  and  administration  elements  of  the 
program  and  for  managing  the  relationships  with  the 
Consortiums  members.  He  manages  the  AHPCRC 
staff  scientists  and  supports  collaboration  activities 
across  the  research  technical  areas,  and  industry/ 
HPC  vendor  partners. 


Barbara  Bryan, 

AHPCRC  Research  and 
Outreach  Manager 

arbara  Bryan,  who  has  managed  AHPCRC 
user  services  at  the  NASA  Ames  facility  since 
it  opened  in  2007,  was  promoted  to  Research 
and  Outreach  Manager  in  2008.  Ms.  Bryan  replaces 
Charles  Peavey,  who  has  assumed  other  duties  at  High 
Performance  Technologies,  Inc.  (HPTi).  Ms.  Bryan 
was  a  staff  member  with  the  first  AHPCRC  Consor¬ 
tium  at  the  University  of  Minnesota  before  coming  to 
HPTi. 


Pritchard,  who  is  based  in  Dayton,  Ohio,  works  for 
High  Performance  Technologies,  Inc.  (HPTi),  the 
administration  and  user  services  member  of  the 
AHPCRC  consortium.  Prior  to  coming  to  HPTi,  he 
was  the  Director  of  Federal  Programs  at  the  Ohio 
Supercomputer  Center.  He  was  previously  employed 
by  Nichols  Research  Corporation  and  BDM 
International. 

Pritchard  is  a  retired  U.S.  Air  Force  Colonel  with  12 
years  of  active  duty  service  and  18  years  as  a  reservist. 
His  technical  background  is  in  computational 
chemistry  and  materials.  He  received  his  B.S.  degree 
in  chemistry  from  Duke  University  and  his  M.S. 
and  Ph.D.  degrees  in  chemistry  from  the  Ohio 
State  University.  He  has  40  years  of  experience  as  a 
researcher,  military  officer,  and  technical/program 
manager  in  defense  research  and  development.  A 


Bryan  works  with  Dr.  Charbel  Farhat,  the  Center 
Director,  and  her  office  is  at  Stanford  University.  She 
works  with  the  Cooperative  Agreement  Manager 
and  Center  Director  to  communicate  research  objec¬ 
tives  and  report  on  accomplishments,  and  align  and 
identify  potential  HPC  resources  on  at-institution  or 
DoD  platforms.  She  also  works  with  each  university’s 
outreach  manager  and  Consortium  representative  to 
establish  an  integrated  outreach  program  emphasizing 
the  Army  research  and  HPC/ computational  science 
objectives. 

Ms.  Bryan  has  extensive  experience  in  high  perfor¬ 
mance  computing  operations  and  support,  parallel 
programming,  scientific  applications,  computer  secu¬ 
rity,  and  scientific  visualization.  She  has  managed  R&D 
projects,  outreach  programs  and  IT  infrastructure  at 
AHPCRC  for  10  years.  A 


Telephone :  (614)  505-6478 
e-mail :  rpritchard@hpti.com 


Telephone  :  (650)  604-3732 
e-mail :  bbryan@hpti.com 
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Cook  Receives  Presidential  Early 
Career  Award  for  Scientists  and 
Engineers 

Ieanine  Cook,  a  principal  investigator  for  the  Army 
High  Performance  Computing  Research  Center 
(AHPCRC),  has  received  the  prestigious  Presiden- 
Early  Career  Award  for  Scientists  and  Engineers 
(PECASE).  President  George  W.  Bush  presented  the 
award  to  Cook  on  December  19,  2008,  at  the  White 
House. 

Cook,  an  associate  professor  at  New  Mexico  State  Uni¬ 
versity’s  Klipsch  School  of  Electrical  and  Computer 
Engineering,  specializes  in  micro-architecture  simula¬ 
tion  techniques,  performance  modeling  and  analysis, 
workload  characterization,  and  micro-architectural 
power  optimizations  for  high  performance  computing 
systems.  She  directs  the  Advanced  Computer  Archi¬ 
tecture  Performance  and  Simulation  Laboratory  at 
NMSU. 

The  PECASE  award  is  one  of  the  highest  honors  that 
the  United  States  government  bestows  on  outstanding 
scientists  and  engineers  beginning  their  independent 
careers.  About  50  U.S.  researchers  receive  PECASE 
awards  each  year,  and  an  individual  can  receive  only 
one  such  award  during  his  or  her  career. 

Cook’s  work  with  the  AHPCRC  consortium  focuses 
on  time-efficient  performance  modeling  and  analy¬ 
sis  of  Army  HPC  applications,  with  the  intention  of 
reducing  the  amount  of  time  and  money  necessary  to 
configure  computer  systems  and  make  procurement 
decisions.  Such  modeling  and  analysis  can  also  as¬ 
sist  in  identifying  optimal  application-to-architecture 
mappings  and  optimizing  code  performance,  operat¬ 
ing  system  services,  and  hardware  design.  She  has 
worked  with  researchers  at  the  University  of  Texas, 

El  Paso,  to  design  and  implement  Chimera,  an  AHP- 
CRC-sponsored  heterogeneous  computing  cluster  that 
links  commercial  multicore  processors,  GPUs,  FPGAs, 
and  accelerators.  Chimera  is  used  to  test  applications 


for  their  ability  to  perform 
well  on  a  variety  of  pro¬ 
cessors,  and  to  identify 
ways  of  reducing  execu¬ 
tion  time  using  specialized 
resources. 

Cook  was  nominated  for 
the  PECASE  award  by 
her  colleagues  at  Sandia 
National  Laboratories  for 
her  work  in  performance 
analysis.  Cook  built  a 
simulator  to  pinpoint  the 
sources  of  Sandia  application 
performance  problems. 

In  an  NMSU  press  release,  Cook,  who  was  paralyzed 
from  the  waist  down  after  a  1982  automobile  acci¬ 
dent,  stated,  “I  want  [the  award]  to  heighten  attention 
for  people  with  disabilities  and  secondly,  I  want  it  to 
heighten  attention  for  young  women.  I  want  young 
women  to  have  role  models  to  say,  ‘If  she  did  that,  I 
can  do  that  too.’  I  really  want  the  world  to  see  that 
people  with  disabilities  are  people — we  are  people— 
we’re  not  to  be  afraid  of.  We  can  be  treated  just  like 
anybody  else.  We’re  not  stupid,  we’re  not  helpless, 
we’re  just  people.” 

In  the  press  release,  NMSU  interim  president  Waded 
Cruzado  praised  Cook,  the  daughter  of  an  Italian 
father  and  a  Hispanic  mother:  “Jeanine  is  a  prime 
example  of  the  outstanding  faculty  talent  New  Mexico 
State  University  strives  to  recruit  and  retain  to  help 
this  land-grant  university’s  mission  and  commitment 
to  providing  the  highest  quality  of  education  for  its 
students.  The  NMSU  family  is  very  proud  of  Jeanine’s 
accomplishments  and  the  honor  and  recognition  that 
she  has  brought  to  herself  and  this  university.  The  dis¬ 
tinction  of  receiving  the  PECASE  award  demonstrates 
how  women,  minorities  and  persons  with  physical 
challenges  can  overcome  the  many  obstacles  placed 
before  them  to  succeed  and  exceed  in  their  careers.”  A 


Jeanine  Cook 


Photo  courtesy  of  New 
Mexico  State  University 
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Scalable  Networks  Technologies 
Receives  Army  Modeling  and 
Simulation  Award 


Rajive  Bagrodia,  president  of  Scalable  Networks 
Technologies,  was  one  member  of  the  Army’s 
Future  Combat  Systems  (Brigade  Combat 
Team)  Communications  Effects  Server  Modeling  and 
Simulation  Team,  which  was  presented  with  an  Army 
Modeling  and  Simulation  Award  at  the  2008  Interser¬ 
vice/Industry  Training,  Simulation  and  Education 
Conference  (ITSEC)  in  Orlando  last  December. 


The  award  recognizes  the  Communications  Effects 
Server  (CES),  an  extension  of  Scalable  Networks  Tech¬ 
nologies’  QualNet  software.  (AHPCRC  is  supporting 
the  deployment  of  QualNet  at  ARL’s  Mobile  Networks 
Modeling  Institute.)  QualNet  and  CES  are  HPC  com¬ 
patible,  which  was  one  reason  that  this  team  won  the 
award,  according  to  Richard  Meyer,  the  company’s 
Director  of  Engineering 

The  CES  is  a  simulation  of  the  Future  Combat  System 
network  transport  layer  that  represents  radios,  ground 


Rajive  Bagrodia  (left),  Richard  Meyer  (right). 
Scalable  Networks  Technologies 

and  airborne  relays,  and  the  Warfighter  Informa¬ 
tion  Network- Tactical  gateways  to  satellite  systems 
connecting  the  FCS  network.  The  CES  can  currently 
simulate  message  delivery  among  at  least  5000  nodes 
running  in  a  near  real-time  environment,  according  to 
the  award  nomination. 

The  nomination  states  that  this  model  runs  in  a  multi¬ 
processor  environment,  and  is  “one  of  the  largest  mo¬ 
bile  ad  hoc  network  models  to  run  in  real  time,  thus 
providing  the  Army  with  a  simulation/ stimulation 
M&S  capability  unequaled  throughout  the  Depart¬ 
ment  of  Defense  or  the  commercial  communications 
industry.”  A 


Multiscale  Ballistic  Fabric  Nets 

continued  from  page  3 


contact  is  distributed  evenly  across  the  entire  prob¬ 
lem  (approximately  one  million  degrees  of  freedom), 
the  algorithm  scales  well  up  to  about  200  processors. 
Optimum  scalability  requires  at  least  50,000  degrees 
of  freedom  per  processor.  The  majority  of  problems 
expected  to  be  of  interest  to  producers  and  users  of 
ballistic  fabrics  will  involve  a  combination  of  local¬ 
ized  and  distributed  contacts,  so  the  scalability  of  the 
algorithm  should  fall  somewhere  between  these  two 
extremes. 


Ballistic  impact  has  been  simulated  for  a  50-caliber 
(0.5-inch)  cylindrical  projectile  and  a  Zylon  square  10 
inches  on  a  side,  fastened  at  the  corners  to  a  substrate, 
for  comparison  with  experimental  studies  done  at 
Berkeley.  In  both  the  experiments  and  the  simulation, 
fabric  failure  occurred  at  the  corners  as  a  result  of 
impact. 

In  the  region  of  the  ballistic  limit  (the  impact  velocity 
below  which  penetration  does  not  occur),  the  simula¬ 
tion  agreed  closely  with  the  experimental  results.  The 
simulation  predicted  a  ballistic  limit  of  39.6  m/sec, 
which  compares  well  with  the  experimentally  ob¬ 
served  ballistic  limit  of  39.3  m/sec.  A 
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Simulation  of  Fracture  and  Penetration 

continued  from  page  5 

Ongoing  Work 

When  many  calculations  are  going  on  all  at  once,  it  is 
necessary  to  build  elements  into  the  computer  code 
that  keep  the  calculations  on  track  and  to  weed  out  so¬ 
lutions  that  are  mathematically  correct,  but  physically 
meaningless.  Toward  this  end,  the  group  is  investigat¬ 
ing  a  new  approach  to  control  the  high-frequency 
velocity  oscillations  that  appear  with  the  adoption  of 
explicit  time  discretizations  such  as  AVI.  These  oscil¬ 
lations  are  enhanced  in  the  presence  of  shocks,  typical 
of  penetration  problems  or  contact  problems.  This 
is  critically  important  for  modeling  the  amount  and 
spatial  distribution  of  energy  absorbed  by  a  ballistic 
gel  when  a  bullet  penetrates  it. 

Raymond  Rickman  is  working  on  a  rigid  bodies  al¬ 
gorithm  that  will  model  the  contact  between  a  bullet 
and  a  ballistic  gel.  He  has  begun  parallel  runs  with 
the  algorithm  using  Army  computing  clusters,  and 
has  shown  that  the  problem  is  scalable.  Ramsharan 
Rangarajan  has  demonstrated  that  a  new  immersed 
boundary  method  for  solid  mechanics  with  nonhomo- 
geneous  boundary  conditions  converges  quadratically 
(a  measure  of  how  fast  a  calculation  reaches  a  solu¬ 
tion),  as  required  by  many  solid  mechanics  problems. 
This  algorithm  is  intended  to  be  part  of  the  automatic 
remeshing  capability  for  the  penetrating  bullet  prob¬ 
lem. 


Mesh  structure  for  projectile  impact  simulation. 


Simulations  for  rigid  body  impact  on  elastic  bodies 
have  been  performed  using  as  many  as  128  proces¬ 
sors  and  3  million  degrees  of  freedom.  Simulations 
for  penetration  have  been  started,  using  a  model  that 
incorporates  a  small  pre-existing  hole  in  the  elastic 
medium.  As  the  work  progresses  further,  a  penetration 
algorithm  will  be  added  so  that  the  pre-opened  hole 
will  not  be  necessary.  A  self-contact  algorithm  is  also 
in  the  works,  which  will  enable  the  simulation  of  the 
cavity  collapsing  after  the  bullet  passes  through. 

The  AHPCRC  group  plans  to  adapt  their  method  to 
model  the  rupture  of  the  material  by  the  projectile. 
They  will  also  develop  the  first  parallel  version  of  the 
penetration  problem  in  gelatin  that  will  run  scalably 
on  HPC  platforms.  This  year,  they  will  perform  the 
first  parallel-explicit  contact  runs  for  the  bullet.  A 
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AHPCRC  Publications  and  Presentations  April  2008-April  2009 

Project  1-1:  Multiscale  Ballistic  Fabric  Nets 

Multi-scale  Modeling  and  Large-scale  Transient  Simulation  of  Ballistic  Fabric.  D.  Powell,  C.  Farhat, ,  T.  I.  Zoh- 
di.  Oral  presentation  CO-03,  Proceedings  of  the  26th  Army  Science  Conference,  Orlando,  FL,  December  2008. 

Project  1-2:  Simulation  of  Fracture  and  Penetration 

•  A  discontinuous-Galerkin-based  immersed  boundary  method.  Adrian  Lew  and  Gustavo  Buscaglia.  Interna¬ 
tional  Journal  for  Numerical  Methods  in  Engineering,  76  (4),  427-454,  2008.  DOI:  10.1002/nme.2312. 

•  Towards  Long-time  Simulation  of  Soft  Tissue  Simulant  Penetration.  Ramsharan  Rangarajan,  Raymond  Ryck- 
man,  Adrian  Lew.  Poster  PP-04,  Proceedings  of  the  26th  Army  Science  Conference,  Orlando,  FL,  December 
2008. 

Project  1-4:  Flapping  and  Twisting  Aeroelastic  Wings  for  Propulsion 

•  Galerkin  Reduced  Order  Models  for  Compressible  Flow  with  Structural  Interaction.  Matthew  F.  Barone, 
Daniel  J.  Segalman,  Heidi  Ihornquist,  Irina  Kalashnikova.  American  Institute  of  Aeronautics  and  Astronautics 
(AIAA)  Paper  No.  2008-0612,  46th  AIAA  Aerospace  Science  Meeting  and  Exhibit,  Reno,  NV,  January  2008). 
(Available  at  www.stanford.edu/~irinak/reno08.pdf) 

•  Flapping  and  Twisting  Aeroelastic  Wings  for  Propulsion.  Humberto  Evans  (presenter),  Mingjun  Wei,  Poster 
presentation,  AHPCRC  exhibit,  Supercomputing  2008,  Austin,  TX. 

•  Flight  of  a  Rufous  Hummingbird  Robotic  Model-PIV  Measurements.  H.  Bocanegra  Evans,  R.  Chavez 
Alarcon,  P.  Ferreira  de  Sousa,  B.  Tobalske,  J.  J.  Allen.  61st  Meeting  of  the  APS  Division  of  Fluid  Dynamics,  San 
Antonio,  TX,  Nov.  23-25,  2008. 

•  A  Computational  Framework  Based  on  Variational  LES  Method  for  the  Multi-disciplinary  Analysis  of  MAVs 
with  Flapping  Wings.  Charbel  Farhat  and  Ajaykumar  Rajasekharan.  Oral  presentation  CO-05,  Proceedings  of 
the  26th  Army  Science  Conference,  Orlando,  FL,  December  2008. 

•  Hamiltonian  structure  and  dynamics  of  a  neutrally  buoyant  rigid  sphere  interacting  with  thin  vortex  rings.  B. 
N.  Shashikanth,  A.  Sheshmani,  S.  Kelly,  and  M.  Wei.  Journal  of  Mathematical  Fluid  Mechanics  (accepted),  2008. 

•  A  Strong- coupling  Approach  to  Simulate  Flexible  Flapping  Wing.  M.  Wei,  T.  Yang,  and  H.  Zhao.  61st  Meeting 
of  the  APS  Division  of  Fluid  Dynamics,  San  Antonio,  TX,  Nov.  23-25,  2008. 

Project  2-1:  Dispersion  of  BWAs  in  Attack  Zones 

•  A  mass,  energy,  vorticity,  and  potential  enstrophy  conserving  lateral  fluid-land  boundary  scheme  for  the  shal¬ 
low  water  equations.  G.  Ketefian,  M.Z.  Jacobson,  EOS  Transactions,  AGU,  89,  Fall  Meeting  Supplement,  2008. 

•  DNS  and  RANS  Simulation  of  Dispersion  Downstream  of  an  Obstacle.  R.  Rossi,  G.  Iaccarino,  Oral  presenta¬ 
tion  QO-04,  Proceedings  of  the  26th  Army  Science  Conference,  Orlando,  FL,  December  2008. 

Project  2-2:  Micro-  and  Nano-fluidic  Devices  for  Sorting  and  Sensing  BWAs 

•  Brownian  demixing  and  wall  effects  in  sedimenting  suspensions  of  orientable  particles.  B.  Hoffman,  E.S.G. 
Shaqfeh,  Physical  Review  E  78,  055301(R),  2008. 

•  The  effect  of  Brownian  motion  on  the  stability  of  sedimenting  suspensions  of  polarizable  rods  in  an  electric 
field.  B.  Hoffman,  E.S.G.  Shaqfeh,  /.  Fluid  Mech.  (in  press,  11/2008). 
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Project  2-2  (continued): 

•  Hydrodynamic  interactions  in  metal  rod-like  particle  suspensions  due  to  induced  charge  electro-osmosis. 

Klint  A.  Rose,  Brendan  Hoffman,  David  Saintillan,  Eric  S.  G.  Shaqfeh,  Juan  G.  Santiago.  Phys.  Rev.  E  79,  011402 
(2009).  doi:  10.1 103/PhysRevE.79.01 1402 

•  Developing  Novel  Empirical  Force  Fields  for  Molecular  Dynamics  Simulations  of  Phase  Transitions  in  BaTiC>3. 
Jose  Solomon  (presenter),  Eric  Darve.  Poster  presentation,  AHPCRC  exhibit,  Supercomputing  2008,  Austin,  TX. 

Project  2-3:  Protein  Structure  Prediction  for  Virus  Particles 

•  Comparison  of  two  approaches  in  the  parallel  computation  of  secondary  structure  topology  determination. 
Saeed  Al  Haj,  Jing  He.  Poster,  IEEE  Bioinformatics  and  Biomedicine  Conference,  Nov  3-5,  Philadelphia. 

•  An  Effective  Convergence  Independent  Loop  Closure  Method  using  Forward-backward  Cyclic  Coordinate 
Descent.  K.  Al  Nasr,  J.  He.  International  Journal  of  Data  Mining  and  Bioinformatics,  3(4),  2009  (in  press). 

•  An  Effective  Convergence  Independent  Loop  Closure  Method  using  Forward/backward  Cyclic  Coordinate 
Descent.  Kamal  Al  Nasr,  Jing  He.  Biosymposium  at  New  Mexico  State  University 

•  Protein  structure  backbone  prediction  for  helical  skeletons  from  low  resolution  protein  density  map.  Kamal  Al 
Nasr,  Weitao  Sun,  Jing  He.  Poster,  IEEE  Bioinformatics  and  Biomedicine  Conference,  Nov  3-5,  Philadelphia. 

•  Improving  ab  initio  protein  Structure  prediction  using  structure  tensor  calculated  from  electron  density  map  of 
a  protein.  Edward  Berndt,  Jing  He.  Biosymposium  at  New  Mexico  State  University 

•  Native  secondary  structure  assignment  has  near  minimum  contact  energy  among  all  possible  assignments.  Jing 
He.  2008  Computational  Structural  Bioinformatics  Workshop,  Philadelphia,  PA,  Nov  3-5. 

•  An  Approximate  Robotics  Algorithm  to  Assemble  a  Loop  Between  Two  Helices.  J.  He,  K.  Al  Nasr.  Proceedings 
of 2007  IEEE  international  Conference  on  Bioinformatics  and  Biomedicine  Workshops,  pp  74-79,  San  Jose,  CA, 
Nov.  2007. 

•  Reduction  of  false  positive  detection  of  helices  from  low  resolution  protein  density  map.  Erick  Rosas,  Jing  He. 
Poster,  IEEE  Bioinformatics  and  Biomedicine  Conference,  Nov  3-5,  Philadelphia. 

•  Protein  structure  prediction  using  geometrical  constraints.  Jing  He.  NMSU-UTEP  CS-MATH  Workshop,  Las 
Cruces,  NM,  April  2008. 

•  Incorporating  Constraints  from  Low  Resolution  Density  Map  in  Ab  Initio  Structure  Prediction  Using  Rosetta. 
Jing  He.  Computational  Structural  Bioinformatics  Workshop,  San  Jose,  CA,  Nov.  2007. 

•  Characterizing  protein  side  chain  interactions  from  the  Protein  Data  Bank.  Joe  Lau,  Jing  He.  Biosymposium  at 
New  Mexico  State  University. 

•  Reduction  of  the  Secondary  Structure  Topological  Space  through  Direct  Estimation  of  the  Contact  Energy 
Formed  by  the  Secondary  Structures.  W.  Sun,  J.  He.  BMC  Bioinformatics,  Vol.  10,  Supplement  1,  S40,  2009. 

(This  was  also  the  topic  of  a  seminar  given  at  NMSU  in  September  2008. 
http://www.cs.nmsu.edu/CSWS/php/seminarDetail.php?colloqia_id=138) 

•  Parallel  Computing  in  Protein  Structure  Topology  Determination.  W.  Sun,  S.  Al-Haj,  J.  He.  Poster  CP-08,  Pro¬ 
ceedings  of  the  26th  Army  Science  Conference,  Orlando,  FL,  December  2008. 

•  Secondary  structure  topology  in  native  protein  is  near  a  minimum  of  contact  energy.  Weitao  Sun,  Jing  He.  Bio¬ 
symposium  at  New  Mexico  State  University 
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Project  2-4:  Nano-Mechanics  of  Metal  Plasticity  of  Thin  Films  and  Cylinders 

•  Surface  controlled  dislocation  multiplication  in  metal  micro-pillars.  C.R.  Weinberger  and  W.  Cai,  Proceedings  of 
the  National  Academy  of  Sciences,  105,  14304-14307,  2008.  doi:10.1073/pnas.08061 18105 

Project  3-1:  Information  Aggregation  and  Diffusion  in  Networks 

•  Enabling  Data  Interpretation  through  User  Collaboration  in  Sensor  Networks.  Eunjoon  Cho,  Kevin  Wong, 
Branislav  Kusy,  and  Leonidas  Guibas.  Proc.  SenSys  Workshop  on  Applications,  Systems,  and  Algorithms  for  Image 
Sensing,  2008. 

•  The  Identity  management  problem —  a  short  survey.  Leonidas  Guibas.  Proc.  International  Conf.  on  Informa¬ 
tion  Fusion,  2008. 

•  Information  Dissemination  and  Aggregation  under  Mobility.  Branislav  Kusy.  Poster  presentation,  26th  Army 
Science  Conference,  Orlando,  FL,  December  2008. 

•  Image  Webs:  Finding  implicit  structure  in  large  image  collections.  Branislav  Kusy.  Poster  presentation,  26th 
Army  Science  Conference,  Orlando,  FL,  December  2008. 

Project  3-2:  Robust  Wireless  Communications  in  Complex  Environments 

•  On  the  Noisy  Interference  Regime  of  the  MISO  Gaussian  Interference  Channel.  B.  Bandemer,  A.  Sezgin,  and  A. 
Paulraj.  CTW  2008  (no  proceedings),  St.  Croix,  US  Virgin  Islands,  USA  ,  May  11-14,  2008. 

Also  presented  at  Asilomar  CSSC  2008,  Pacific  Grove,  October  26-29,  2008. 

•  Antenna  Selection  in  Space-Time  Block  Coded  Systems:  Performance  Analysis  and  Low-Complexity  Algo¬ 
rithm.  C.-Y.  Chen,  A.  Sezgin,  J.M.  Cioffi,  and  A.  Paulraj.  IEEE  Transactions  on  Signal  Processing,  56(7),  3303- 
3314,  2008.  doi:  10.1109/TSP.2008.917856 

•  Passive  Sensor  Imaging  Using  Cross  Correlations  of  Noisy  Signals  in  a  Scattering  Medium.  Josselin  Gamier  and 
George  Papanicolaou.  Eos  Transactions,  89(53),  American  Geophysical  Union  2008  Fall  Meeting  Supplement, 
Abstract  S23D-04. 

•  Feedback  Reduction  in  Uplink  MMO  OFDM  Systems  by  Chunk  Optimization.  E.A.  Jorswieck,  A.  Sezgin,  B. 
Ottersten,  and  A.  Paulraj.  EURASIP  Journal  on  Advances  in  Signal  Processing,  Special  Issue  °MIMO  Transmission 
with  Limited  Feedback”  ,  Vol.  2008,  Article  ID  597072.  doi:10.1 155/2008/597072 

•  Feedback  Reduction  in  Uplink  MIMO  OFDM  Systems  by  Chunk  Optimization.  E.A.  Jorswieck,  A.Sezgin,  B. 
Ottersten,  and  A.  Paulraj.  ICC  2008,  Beijing,  China,  May  19-23,  2008. 

•  Guaranteed  Performance  Region  in  Fading  Orthogonal  Space-Time  Coded  Broadcast  Channels.  E.A.  Jor¬ 
swieck,  A.  Sezgin,  B.  Ottersten,  and  A.  Paulraj.  EURASIP  Journal  on  Wireless  Comm,  and  Networking,  Special 
Issue:  Theory  and  Applications  in  Multiuser/Multiterminal  Communications,  Vol.  2008,  Article  ID  268979. 
doi:10.1 155/2008/268979 

•  Framework  for  analysis  of  opportunistic  schedulers:  average  sum  rate  vs.  average  fairness.  E.  Jorswieck, 
A.Sezgin,  andX.  Zhang.  RAWNET  2008,  Berlin,  Germany,  March  31,  2008.  (http://www.ee.kth.se/php/modules/ 
publications/ reports/2008/IR-EE-SB_2008_0 11. pdf) 

•  Opportunistic  Scheduling  for  Multiantenna  Cellular:  Interference  Limited  Regime.  S.  Pereira,  A.  Paulraj  and 
G.  Papanicolaou.  Signals,  Systems  and  Computers,  2007.  ACSSC  2007.  Conference  Record  of  the  41st  Asilomar 
Conf.  on  Signals,  Systems,  and  Computers,  Pacific  Grove,  CA,  November  4-7,  2007.  pp  1171-1175,  doi  10.1109/ 
ACSSC.2007.4487408 

•  Interference  Limited  Broadcast:  Role  of  Interferer  Geometry.  S.  Pereira,  A.  Sezgin,  A.  Paulraj,  and  G.  C.  Papani¬ 
colaou.  Information  Theory,  2008.  ISIT  2008.  IEEE  International  Symposium,  Toronto,  Canada,  July  6-11,  2008. 
pp  757-761,  doi  10.1109/ISIT.2008.4595088 
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Project  3-2  (continued): 

•  Interaction  between  scheduling  and  user  locations  in  an  OSTBC  coded  downlink  system.  A.  Sezgin,  E.A.  Jor- 
swieck,  M.  Charafeddine,  A.  Paulraj.  7th  International  ITG  Conference  on  Source  and  Channel  Coding,  Ulm, 
Germany,  January  2008. 

•  On  the  relation  of  OSTBC  and  code  rate  one  QSTBC:  Average  Rate,  BER  and  Coding  Gain.  A.  Sezgin,  E.A. 
Jorswieck,  O.  Henkel,  S.  Pereira,  and  A.  Paulraj.  IEEE  Transactions  on  Signal  Processing,  56(10),  4879-4891, 
2008.  doi:  10.1109/TSP.2008.923202 

•  Tile  based  MIMO  OFDMA  systems:  Impact  of  outdated  feedback.  Sezgin,  E.A.  Jorswieck,  B.  Bandemer,  and 
A.  Paulraj.  Presented  at  Asilomar  CSSC  2008,  October  26-29,  2008. 

•  Multicell  Optimization  for  Diversity  and  Interference  Mitigation.  O.  Oteri,  A.  Paulraj.  IEEE  Transactions  on 
Signal  Processing,  56,  2050-2061,  2008. 

•  Where  to  place  interferers  in  a  wireless  network.  A.  Sezgin,  E.  Jorswieck,  and  A.  Paulraj.  Information  Theory 
Workshop,  2008.  ITW  ‘08.  IEEE,  Porto,  Portugal,  May  5-9,  2008.  pp  348-352,  doi  10.1109/ITW.2008.4578683 

•  Complete  Characterization  of  the  Equivalent  MIMO  Channel  for  Quasi-Orthogonal  Space-Time  Codes.  A. 
Sezgin  and  T.J.  Oechtering.  IEEE  Transactions  on  Information  Theory,  54(7),  3315-3327,  2008.  doi:  10.1109/ 
TIT.2008.924705 

•  Mobility  Dependent  Feedback  Scheme  for  point-to-point  MIMO  Systems.  G.  Vazquez- Vilar,  V.  Majjigi,  A. 
Sezgin,  and  A.  Paulraj.  Presented  at  Asilomar  CSSC  2008,  October  26-29,  2008. 

•  Statistical  Opportunistic  Scheduling  With  Tap  Correlation  Information  for  an  OFDMA  System  in  Uplink.  E. 
Yoon,  D.  Tujkovic,  A.  Paulraj.  IEEE  Transactions  on  Signal  Processing,  57, 1708-1714,  2008. 

•  Patent  disclosure  for  “Signal  Archiving  in  Ad  hoc  Networks”  made  to  Stanford  University,  lists  AHPCRC  as 
funding  source,  thus  giving  certain  rights  for  the  U.S.  Government. 

Project  3-3:  Secure  Sensor  Data  Dissemination  and  Aggregation 

Impact  of  Distributed  Jamming  on  Sensor  Network.  Hong  Huang,  Amjad  Abu-Baker.  Poster  presentation,  26th 
Army  Science  Conference,  Orlando,  FL,  December  2008. 

Project  4-3:  Simulation  &  Modeling  to  Enhance  the  Performance  of  Systems  of  Multicore  Processors 

•  Analysis  of  a  Checkpoint/Restart  Model  to  Uncover  Properties  that  Can  Lead  to  Improved  Application  Per¬ 
formance.  S.  Arunagiri.  The  4th  Joint  NMSU/UTEP  Workshop  on  Mathematics,  Computer  Science,  and  Com¬ 
putational  Sciences,  University  of  Texas  at  El  Paso,  November  8,  2008. 

•  Opportunistic  Checkpoint  Intervals  to  Improve  System  Performance.  Sarala  Arunagiri,  John  T.  Daly,  Patricia 
J.  Teller,  Seetharami  Seelam,  Ron  A.  Oldfield,  Maria  Ruiz  Varela,  Rolf  Riesen.  Technical  Report  UTEP-CS-08-24, 
April  2008.  http://www.cs.utep.edU/vladik/2008/abstr08.html#24 

•  Monte  Carlo  Processor  Modeling.  J.  Cook.  Los  Alamos  Computer  Science  Symposium,  Workshop  on  Perfor¬ 
mance  of  Extreme-Scale  Systems  and  Applications,  October  13-15,  2008,  Santa  Fe,  NM. 

•  Modeling  Processor  Performance  Using  Monte  Carlo  Simulation.  J.  Cook,  D.  Balise.  Poster  presentation, 
AHPCRC  exhibit.  Supercomputing  2008,  Austin,  TX. 

•  Improving  Throughput  of  SMT  Processors  using  Application  Signatures  and  Thread  Priorities.  Mitesh 
Meswani.  Supercomputing  2008,  December  1-4,  Austin,  Texas. 
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•  Measuring  and  Validating  Metrics  used  to  Estimate  Microarchitecture  Resource  Utilization:  A  Case  Study 
of  the  IBM  POWER5  Processor.  M.  R.  Meswani,  P.  J.  Teller,  S.  Arunagiri.  Proceedings  of  the  2008  Live-Virtual- 
Constructive  Conference  (hosted  by  ITEA  White  Sands  Chapter),  January  12-15,  2009,  El  Paso,  TX.  Accepted  for 
publication. 

•  A  Study  of  the  Influence  of  the  POWER5  Dynamic  Resource  Balancing  Hardware  on  Optimal  Hardware 
Thread  Priorities.  P.  C.  Trillo,  M.  R.  Meswani,  P.  J.  Teller,  S.  Arunagiri.  Proceedings  of  the  2008  Live-Virtual- 
Constructive  Conference  (hosted  by  ITEA  White  Sands  Chapter),  January  12-15,  2009,  El  Paso,  TX.  Accepted  for 
publication.  Selected  as  the  Best  Undergraduate  Student  Paper. 

Project  4-4:  High  Performance  Optimization  Library 

•  A  Large-Scale  Quadratic  Programming  Solver  Based  on  Block- Lu  Updates  of  the  KKT  System.  Hanh  M. 
Huynh.  Doctoral  dissertation,  September  2008.  (available  at  http://www.stanford.edu/group/SOL/dissertations/ 
hanh-thesis.pdf) 

•  An  algorithm  for  solving  nonlinear  optimization  problems  with  linear  constraints  and  binary  variables.  Kien 
Ming  Ng,  Walter  Murray.  IMA  workshop  on  Mixed-Integer  Nonlinear  Optimization:  Algorithmic  Advances  and 
Applications,  November  17-21,  2008,  University  of  Minnesota.  (http://www.ima.umn.edu/2008-2009/SWll.17- 
2 1 .08/abstracts.html#7322) 

Also  accepted  for  publication  by  the  Journal  of  Computational  Optimization  and  Applications. 

Project  4-5:  PFMMPACK 

Fourier  Based  Fast  Multipole  Method.  Cris  Cecka,  Eric  Darve.  Poster  QP-04,  Proceedings  of  the  26th  Army  Sci¬ 
ence  Conference,  Orlando,  FL,  December  2008. 

Project  4-6:  Hybrid  Optimization  Schemes  for  Parameter  Estimation  Problems 

•  Numerical  Testing  of  Parameterization  Schemes  for  Solving  Parameter  Estimation  Problems.  M.  Argaez,  L. 
Velazquez.  Poster  CP-06,  Proceedings  of  the  26th  Army  Science  Conference,  Orlando,  FL,  December  2008. 

•  Compressive  Sensing  Reconstruction  with  Prior  Information  by  Iteratively  Reweighted  Least  Squares.  M. 
Argaez,  L.  Velazquez,  C.  Quintero,  R.  Von  Borries,  C.  Potes,  C.  Mendes.  IEEE  Transactions  on  Signal  Processing 
(in  press). 

•  A  Sequential  Quadratic  Program  for  Nonconvex  Lp-Regularized  Metric  Problems.  Miguel  Argaez.  Interna¬ 
tional  SCAN  (Scientific  Computing,  Computer  Arithmetic,  and  Verified  Numerical  Computations)  Conference, 
UTEP,  El  Paso,  TX,  September  2008. 

•  Poster  presentation  (title  not  available).  Carlos  Quintero.  NSF  CoPROD  (Constraint  Programming  and  Deci¬ 
sion  Making)  Workshop,  UTEP,  El  Paso,  TX,  October  2008.  (http://coprod.cs.utep.edu) 

•  Parallel  Global  Optimization  Schemes  for  Solving  Parameter  Estimation  Problems.  Carlos  Quintero.  Poster 
presentation,  AHPCRC  exhibit,  Supercomputing  2008,  Austin,  TX. 

•  Parameter  optimization:  Algorithms  combine  multiple  computing  methods,  computer  architecture  innova¬ 
tions.  Carlos  Ramirez.  Poster  presentation,  AHPCRC  exhibit.  Supercomputing  2008,  Austin,  TX. 

•  Poster  presentation  (title  not  available).  Carlos  Quintero.  SACNAS-UTEP  Expo  2008,  El  Paso,  TX,  May  2008. 

•  Hybrid  methods  for  solving  large-scale  parameter  estimation  problems.  Leticia  Velazquez,  Carlos  Quintero 
(presenter),  and  Miguel  Argaez.  International  SCAN  Conference,  UTEP,  El  Paso,  TX,  September  2008.  (http:// 
www.cs.utep.edu/interval-comp/scan08quintero.ppt) 

•  Invited  talk  (title  not  available).  L.  Velazquez.  NSF  CoPROD  Workshop,  UTEP,  El  Paso,  TX,  September  2008. 
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